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Torque Scaling in Turbulent Taylor-Couette Flow with Co- and Counterrotatinj 

Cylinders 
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We analyze the global transport properties of turbulent Taylor-Couette flow in the strongly 
turbulent regime for independently rotating outer and inner cylinders, reaching Reynolds numbers 
of the inner and outer cylinders of Re; = 2 x 10 6 and Re = ±1.4 x 10 6 , respectively. For all Re;, 
Re , the dimensionless torque G scales as a function of the Taylor number Ta (which is proportional 
to the square of the difference between the angular velocities of the inner and outer cylinders) 
with a universal effective scaling law G oc Ta 0,88 , corresponding to Nu^ oc Ta ' 38 for the Nusselt 
number characterizing the angular velocity transport between the inner to the outer cylinders. 
The exponent 0.38 corresponds to the ultimate regime scaling for the analogous Rayleigh-Benard 
system. The transport is most efficient for the counterrotating case along the diagonal in phase 
space with uj » — 0.4u;;. 
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Global transport properties of turbulent flows are of 
prime importance for many applications of fluid dynam- 
ics, but also for a fundamental understanding, as they 
reflect the interplay between the boundary layer and 
the bulk. The two canonical systems used to analyze 
the transport properties in closed turbulent systems are 
Rayleigh-Benard (RB) convection and Taylor-Couette 
(TC) flow, and they are conceptually closely related 0- 
3|. In RB flow, heat (in dimensionless form the Nusselt 
number) is transported from the hot bottom plate to the 
cold top plate [H, j whereas in TC flow angular velocity 
is transported from the inner to the outer cylinder (for 
LJi > uj ). In analogy to RB flow, Eckhardt et al. Q 
identified, from the underlying Navier-Stokes equations, 



r 3 ((u r uj} A t - vd r (u) A<t ), 



(1) 



as a relevant conserved transport quantity, representing 
the flux of angular velocity from the inner to the outer 
cylinder. Here u r (u^,) is the radial (azimuthal) velocity, 
uj = u^/r the angular velocity, and {. . .) At characterizes 
averaging over time and an area with constant r from 
the axis. J" is made dimensionless with its value Jf am = 
2vrfr1(uji — uj )/(r 2 — rf) for the laminar case, giving a 
"Nusselt number" as a dimensionless transport quantity, 



n Uu = rjj? 



i am ' 



(2) 



where r; i0 and w; i0 denote the radius and the angular 
velocity of the inner and outer cylinders, respectively, and 
v is the kinematic viscosity of the fluid. Nu w is closely 
connected to the torque r that is necessary to keep the 
inner cylinder rotating at a constant angular velocity or, 
in dimensionless form, to 



G 



Nil 



•J i. 



where I is the height of the cylinder and pfl u id the den- 
sity of the fluid. Yet another often used possibility 
to represent the data is the friction coefficient c/ = 
[(1 - ?? ) 2 A]G/Re? |. 

For RB flow, the scaling properties of the Nusselt num- 
ber in the fully turbulent regime (i.e., for very large 
Rayleigh numbers, say Ra > 10 10 ) have received tremen- 
dous attention in the last decade and various heat flux 
measurements have been performed; again, see the review 
article Q. In contrast, TC flow in the fully turbulent 
regime has received much less attention, with the only 
exception being the Texas experiment by Swinney, Lath- 
rop, and coworkers yM^. In that experiment a Reynolds 
number Re.; = 10 6 of the inner cylinder was reached (with 
the outer cylinder at rest) and an effective power law of 
G oc Re] with 7 w 1.6 - 1.86 was detected @, Q in 
the turbulent regime, though the scaling properties are 
not particularly good and, strictly speaking, 7 depends 
on Re^; i.e., there is no pure scaling. Indeed, in Refs. 



lOj, lll| we have argued that there should be a smooth 



Nu w G, 



lam 1 



(3) 



transition from 7 = 3/2 for the small Reynolds number 
of a boundary layer dominated flow to 7 = 2 for the 
larger Reynolds number of a flow dominated by the tur- 
bulent bulk. Turbulent TC experiments get even more 
scarce for TC flow with inner and outer cylinders ro- 
tating independently. We are only aware of the Wendt 
experiments in the 1930s [l2| . reaching Re; « 10 5 and 
Re D w ±10 5 and the recent ones by Ravelet et al. [3], 
reaching Re; «5x 10 4 and Re — ±2x 10 4 . The hitherto 
explored phase diagram of TC flow with independently 
rotating cylinders is shown in Fig. HJa). An alternative 
representation of the phase diagram is given by Dubrulle 
et al. 15], who introduce a shear Reynolds number and 
a rotation number as an alternative representation of the 
phase space (see below). 
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FIG. 1: (a) Explored phase space (Re , Rei) of TC flow with independently rotating inner and outer cylinders. To the right 
of the horizontal axis the cylinders are corotating, to the left of it they are counterrotating, and a log-log representation has 
been chosen. Experimental data by Wendt ^| (pluses), Taylor [2lJ (left triangles), Smith and Townsend [221 (open circles), 
Andereck et al. 23] (grey box), Tong et al. 13| (upward triangles), Lathrop et al. [g] (stars), Ravelet et al. Jpjj] (hexagrams), 
Borrero-Echeverry et al. [24} (upward solid triangles), and simulations by Pirro and Quadrio [2^| (solid squares), Bilson and 
Bremhorst 26] (downward solid triangles), and Dong [27l. l2Sl] (solid circles). The dashed lines are Esser and Grossmann's [2S| 
estimate for the onset of turbulence with n = 0.71. The many data points in the small Reynolds number regime of pattern 
formation and spatial temporal chaos (see e.g. 23, 30, 31]) have not been included in this phase diagram. Our data points for 
this publication are the black diamonds. Our data points for this publication are the black diamonds, (b) Our data points in 
the phase diagram on a linear scale, (c) Our data points in the phase space (Ta, a); note that Ta also depends on a. (d) Our 
data points in the phase space (Ta, Rn) [see Eq. J?}]. 



In the phase diagram [Fig.QJa)] we have also added the 
Rej, Re Q numbers that we explored with our newly con- 
structed Twente turbulent TC facility (T 3 C), which we 
have described in great detail in Ref. [lfj. In short, at this 
facility the cylinders are £ — 0.97 m high; the inner cylin- 
der has a radius of = 20 cm, the outer cylinder has a 
radius of r Q — 27.9 cm, and the maximal inner and outer 
angular velocity are Wi/27r = 20 Hz and uj /2it = ±10 Hz, 
respectively, corresponding to Re^ — riUJidjv = 2 x 10 6 
and Re Q = r uj d/v = ±1.4 x 10 6 , with d = r Q - ri. The 
system is fully temperature controlled through cooling of 
the upper and lower plates. The torque is measured at 



the middle part of the inner cylinder (similar to [7|) by 
load cells imbedded inside the inner cylinder and not by 
measuring the torque through the seals. One of the goals 
we want to achieve with this new facility is to explore 
the (Rei, Re , r\ = ri/r ) parameter space, thus entering 
terra incognita, and measure the torque [i.e., transport 
of the angular velocity (Nu w in dimensionless form) or, 
again expressed differently, the overall drag] and the in- 
ternal Reynolds number of the flow. 

In this Letter we will focus on the required torque for 
fully developed turbulent flow (Re^, Re„ > 10 5 ), where 
7/ = 0.716 with independently rotating inner and outer 
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cylinders, which hitherto has not been explored. The 
examined parameter space in this Letter is shown in the 
space of (Rej, Re D ) in Fig. 1(b), (Ta, a) in Fig. 1(c), and 
(Ta, Rq) in Fig. 1(d), to be explained below. We will 
not address the question whether pure scaling laws exist: 
First, the explored Reynolds number range is too short to 
answer this question, and second, the earlier work 0,0-11, 
IqL 11 1 gives overwhelming experimental and theoretical 
evidence that there are no pure scaling laws, even up to 
Reynolds numbers of 10 6 . So all scaling exponents in this 
Letter have to be read as effective scaling laws. 

Our results for the counterrotating case for the dimen- 
sionless torque G as function of Re^ for fixed Re are 
shown in Fig. [2] One immediately sees that counterro- 
tation enhances the torque (and thus the overall drag), 
but that for general Re Q ^ the effective power law 
G oc Re*' 76 , that holds in the case of inner cylinder ro- 
tating only, gets lost; in fact, there is even no effective 
power law at all. 

How can we represent the data to better reveal the 
transport properties of the system? The analysis of Eck- 
hardt et al. [3J and the analogy of the TC system to the 
RB system suggest to better plot Nu^ as a function of 
the Taylor number 



Ta = -ad (n + r Q ) (u>i - uj q ) v 



2, ,-2 



(4) 



where a = {[(1 + ij) / '2] / \/rj} 4 , i.e., along the diagonals 



(jj n = — acj; 



(5) 

in the parameter space [17|, Fig. QJb). Indeed, Eck- 
hardt et al. derived, from the underlying Navier- 
Stokes equation, the exact relation 



^ 3 d- 4 o-- 2 Ta(Nu w - 1) 



(6) 
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FIG. 2: The dimensionless torque G(Rei) for counterrotating 
TC flow for four different fixed values of Re = —1.4 x 10 6 , 
—0.8 x 10 6 , —0.4 x 10 6 , and (top to bottom data sets); see 
inset for the probed area of the parameter space. 



for the excess kinetic energy dissipation rate e w [i.e., 
the total kinetic energy dissipation rate e minus the ki- 
netic energy dissipation rate in the laminar case e; am = 
Avr 2 r 2 (r i +r )~ 2 d~ 2 (u l —uj ) 2 ]. In Eq. ([6} a can be inter- 
preted as a (geometric) Prandtl number, and Ta and Nu u 
are the exact TC analogs to the Rayleigh and Nusselt 
numbers in RB flow. Along the diagonal, Eq. (|5J), in pa- 
rameter space, one has Ta = \ad 2 {ri + r ) 2 (l + a) 2 w 2 v~ 2 , 
and the well-studied [H effective scaling law Nu ex Ra 7 
for RB flow (with 7 0.31 [4, 18]) would now correspond 
to an effective scaling law Nu w oc Ta 7 for TC flow. 

Nu w vs Ta is shown in Fig.^a) for various a, i.e., along 
various straight lines through the origin of the parameter 
space, Fig. QJb). A universal, i.e. a-indepedent, effective 
scaling Nu w cx Ta 7 with 7 ps 0.38 is clearly revealed. 
This corresponds to a scaling of G oc Re^' 76 for the di- 
mensionless torque along the straight lines, Eq. (JS), in 
the parameter space, Fig.Q] to c/ oc Re~ ' 24 for the drag 
coefficient, and to G oc Ta ' 88 . The compensated plots 
Nu^/Ta 38 in Fig. G^b) demonstrate the quality of the 
effective scaling and, in addition, show the a dependence 
of the prefactor of the scaling law. 



10" 



+ 


a = 2 


> 


0.14 




1.00 


O 


0.00 




0.71 


□ 


-0.14 




0.40 





-0.20 




0.20 


V 


-0.40 




FIG. 3: (a) Nu u vs Ta for various a; see Fig. []Jb) for the 
location of the data in parameter space. A universal effective 
scaling Nu^ oc Ta ' 38 is revealed. The compensated plots 
Nu^/Ta ' 38 in (b) show the quality of the effective scaling 
and the a-dependent prefactor of the scaling law. 
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FIG. 4: Prefactor of the effective scaling law Nu„ oc Ta 
(shown in Fig. \3§ as a function of a — —uj /uji. The inset 
shows the effective exponents 7 which results from an indi- 
vidual fit of the scaling law Nu„ oc Ta 7 . 



The a dependence of the prefactor Nu w /Ta ' 38 is plot- 
ted in Fig. 21 It shows a pronounced maximum around 
a = 0.4, i.e., for the moderately counterrotating case, sig- 
naling the most efficient angular velocity transport from 
the inner to the outer cylinder at that value. We men- 
tion that it is obvious that this curve has a maximum, 
as in both limiting cases a ±00 (rotating of the outer 
cylinder only) the flow is laminar and Nu w = 1, but it is 
interesting to note that the maximum does not occur for 
the most pronounced counterrotating case lo = ~Ui (or 
a = 1). Compared to the case of pure inner cylinder ro- 
tation (a = 0), at a = 0.4 the angular velocity transport 
from the inner to the outer cylinder is enhanced by more 
than 20%. 

The parameter a = —ui /u>i is connected to the so- 
called rotation number 



Rq = (1 - 77) (Re, + Re )/(7yRe - R Ci 



(7) 



introduced by Dubrulle et al. [15] and used by Ravelet et 
al. [LJ], namely, Rn = {l—rj)r]^ 1 {a—rj){a+l)~ 1 . We also 
plot our data points in the phase space of (Ta, Rn), as 
shown in Fig. 1(d). The optimal value a ps 0.4 we found 
for the transport properties of the system corresponds 
to Rq s» —0.09. In this Letter we prefer a as compared 
to as the sign of a immediately signals whether the 
system is corotating or counterrotating. 

In conclusion, we have explored the terra incognita of 
fully developed turbulent TC flow with independently ro- 
tating inner and outer cylinders, beyond Reynolds num- 
bers of 10 6 , finding a universal scaling law G oc Ta ' 88 , 
corresponding to Nu w 



Ta 1 '" 18 , for all (fixed) a = 
—LUo/uii, with optimal transport quantities at a « 0.4. It 
is remarkable that the effective scaling exponent 0.38 ex- 
actly resembles the analogous effective scaling exponent 
in Nu oc Ra ' 38 in RB convection in the ultimate regime 



of thermal convection 19|, 2(| , reflecting the analogy be- 
tween TC and RB flow also in the strongly turbulent 
regime. 

The next steps will be to further extend the parame- 
ter space, Fig. [TJ towards further radius ratios n to see 
whether the observed universality carries on towards an 
even larger parameter range, and to also measure the 
Taylor-Reynolds number and the wind Reynolds numbers 
of the internal flow, which are closely connected to Nu^ 
and for which theoretical predictions exist [3| . With such 
measurements and characterizations of the flow struc- 
tures, we will also be able to check whether these are 
reflected in the overall transport properties. 
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ous stimulating discussions about TC flow over the years, 
D. Narezo Guzman for help with some measurements, 
TCO-TNW Twente for the collaboration on building up 
the facility, and STW, which is financially supported by 
NWO, for financial support. 

Note added in proof. — After submission of our Letter, 
Dan Lathrop mase us aware of the parallel work |32j |. 
independently confirming the peak in the dimensionless 
torque as a function of a. 
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